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Abstract. In this paper we develop a technique for determining interatomic 
potentials in materials in the quantum regime from single-shell Extended 
X-ray Absorption Spectroscopy (EXAFS) spectra. We introduce a 
pair distribution function, based on ordinary quantum time-independent 
perturbation theory. In the proposed scheme, the model potential parameters 
enter the distribution through a fourth-order Taylor expansion of the potential, 
and are directly refined in the fit of the model signal to the experimental 
spectrum. We discuss in general the validity of our theoretical framework, 
namely the quantum regime and perturbative treatment, and work out a 
simple tool for monitoring the sensitivity of our theory in determining lattice 
anharmonicities based on the statistical F-test. As an example, we apply our 
formalism to an EXAFS spectrum at the Ag K-edge of Agl at T = 77 K. 
We determine the Ag-I potential parameters and find good agreement with 
previous studies. 



1. Introduction 

It is well known that EXAFS is a very sensitive and accurate technique for probing 
the structural and dynamical properties of materials in the neighborhood of the 
photoabsorber atom [1]. In particular, the damping of the EXAFS signal induced 
by thermal broadening in the distribution of absorber-neighbour distances carrys 
quantitative information on the corresponding interatomic potentials [2]. 

The effects of thermal disorder are usually accounted for by introducing 
a temperature-dependent pair distribution function (PDF) g(r,T), such that 
g s (r,T)dr is the normalised probability that the distance of an atom in the the 
s-th shell from the absorber lies in the interval [r,r + dr] at the temperature T. The 
corresponding single-shell EXAFS Xs(k) is then obtained as 



It is clear that the absorber-neighbour effective potential determine the shape of the 
PDF. The simplest structural model one can introduce is represented by a Gaussian 
distribution, i.e. harmonic interatomic potential. In this case, integration of eq. (1) 
is straightforward and the result is the well-known Debye-Waller multiplicative 
damping factor e~ 2a k , where a 2 = ([r — (r)} 2 ) = g(r,T)[r — (r)] 2 dr. 

The higher-order terms in the potential expansion are usually accounted for 
via the cumulant method [3] . The cumulants C n are defined as the coefficients that 
enter the MacLaurin expansion of the function ln[F(2fc)], where F(2k) is the Fourier 
transform of the effective distribution function G(r, A) = g(r) e~ 2r ^ x /r 2 , A being 
the photoelectron mean free path. The utility of this method is that the cumulants 
are related to the moments of the effective distribution. In particular, the first 
cumulant C\ is the mean value of the interatomic distance, while C2 is the variance 
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of the effective distribution. The higher-order cumulants C 3 and C4 are related 
to the skewness (asimmetry) and kurtosis (deviation from the Gaussian shape) of 
the distribution, respectively. The cumulant-expansion technique has been applied 
succesfully in the case of moderate anharmonicities both to bulk materials [4] and 
to the study of surfaces [5, 6]. 

In such a framework, the following step towards a detailed understanding of 
disordered systems is to establish the direct relationship between the interatomic 
potential and cumulants. In the case of harmonic crystals, the second-order 
cumulants have been calculated quantum-mechanically for Debye crystals [7] and 
simple molecules [8]. Other attempts in this direction exist in the literature which 
include third-order anharmonicities for simple systems [9, 10]. 

A more straightforward alternative approach is to directly integrate the EXAFS 
by calculating explicitely the whole pair distribution function corresponding to a 
certain model potential. The simplest way a pair potential V(r) may directly enter 
the distribution function is through the classical configurational integral scheme 

e -V(r)/k B T 

9^T)= x . (2) 

/ e -V(r)/k B T dr 

Jo 

where k B is the Boltzmann constant. Expression (2) has been applied to the study of 
metals and ionic systems with different model potentials, from Lennard-Jones [11], 
Morse [12] to generic three-parameter Taylor expansions [13]. However, eq. (2) is 
based on a classical treatment of the atomic vibrations. In general, depending both 
on the temperature and potential stiffness, the classical approximation may break 
down. In this case the full quantum treatment of lattice dynamics is in order, and 
one can proceed as follows. 

First the Schrodinger equation of the absorber-neighbour pair has to be solved, 
and the eigenvectors tjj n and eigenvalues E n computed. The quantum-mechanical 
analogue of expression (2) can then be written as 

n 

where r = xi — X2I is the spatial coordinate describing the relative motion of the 
pair. Eq. (3) has been first introduced and used in ref. [14] to study the interatomic 
potential of the Cu-0(4) pair in the YBCO superconductor. The procedure is 
the following: one first introduces a model potential which is characterised by 
a set of parameters {A}. Then a numerical routine is set up, which solves the 
radial Schrodinger equation for a particular choice of the set {A}, builds the radial 
distribution function (3) and calculates the EXAFS signal x(fc, {A}). This routine 
is then incorporated in the fitting program that refines the free parameters {A} on 
a set of experimental data through ordinary x 2 minimisation. 

This procedure has the advantage that it allows an arbitrary analytical 
potential function to be used. This is the case of ref. [14], where a double-well 
potential is found. However, it is rather cumbersome and of little practical utility 
for routine fittings. In particular, a more concise and handy way of calculating 
expression (3) in a closed form would be of great advantage. 

In this paper we compute an analytical expression for the function (3), based 
on a simple Taylor expansion of the potential function, and cast it in a simple form, 
suitable for inclusion in a simple routine attached to the fitting machinery. We 
organise our paper as follows. In section 2 we develop our quantum pair distribution 
function (QPDF). In section 3 we discuss on general grounds the validity of our 
theoretical framework. Furthermore, we develop a statistical tool to assess the 
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sensitivity of the QPDF to the parameters describing the anharmonicities in the 
potential. Finally, in section 4 we test the QPDF by fitting an EXAFS spectrum of 
Agl at T = 77 K. We end the paper by summarising our results and drawing our 
conclusions in section 5. 



2. The Quantum Pair Distribution Function 



Let us consider the pair formed by the photoabsorber and one of its neighbours 
from a given coordination shell. Let Xi and (i = 1, 2) denote the position vectors 
and atomic masses of the two atoms, respectively. Let V(r) be the corresponding 
interatomic potential. We can write its Taylor expansion in the following fashion 



V(r) = -k 2 (r - r ) 2 + V^r) + 0(\r - r \ 5 ) 



(4) 



where 



Vi(r) = -h(r- 



-r ) 3 + -fc 4 (r-r ) 4 



and 



(m-1)! 



d m V(r) 



dr r 



,(5) 

r=r 

r being the equilibrium interparticle distance, given by the condition 
[dV {r) / dr] r=ro = 0. We can follow the ordinary procedure to first separate the 
two-body Schrodinger equation by introducing the relative and centre-of-mass 
coordinates, and then decouple the angular and radial degrees of freedom in the 
Schrodinger equation for the radial motion. The pair wavefunction then reads 



*(xi,x 2 ) = Vg(X) 



u{r) 



(6) 



where X = (mixi +J7i2X2)/(mi -\-rn2), r = |xi — X2 and Y" ; m (#, (f>) are the spherical 
harmonic functions [15]. We require the wavefunction of the pair to have spherical 
symmetry, since we do not want the PDF to depend on the orientation of the 
absorber-neighbour bond in the crystal \. We therefore set I = 0. The radial 
equation then reduces to the one-dimensional problem 



H 2 d 2 u 



+ 



h 2 (r 



rof + V^r) 



Eu 



(7) 



2fi dr 2 

where \i is the reduced mass of the pair and we require u(r)\ r=0 = 0. In the 
spirit of ordinary time-independent perturbation theory, we consider the harmonic 
Hamiltonian as the unperturbed problem and the potential V\ as the perturbation. 

It is convenient to adopt the formalism of second quantisation. The 
unperturbed problem is defined by the eigenvectors \n) and the corresponding 
eigenvalues En" 1 — Tiuj[n + 1/2] (n = 0,1,2,...), where lo = \Jk2l y,- Recalling 
the well known commutation relations between creation and annihilation operators 
d^ and a, respectively, it is straightforward to write down the expression of the 
perturbation potential (5) in the n-representation. We get 



V = Tiuj A3 [a 



where 



t3 



+ lind ] + 3(n + l)a] + 
hujK i [o) A + a 4 + 2(2n - l)a t2 + 2(2n + 3)a 2 + 3(n + l) 2 + 3n 2 



(8) 



A m = 



1 



, . ( (m = 3,4) and x = ^n/Jiu) .(9) 

2 m ' z m \ Tllu 

Recalling the definition of creation and annihilation operators we can easily evaluate 
the matrix elements (fc|V|n). Let us define the four-dimensional quantity /3 as 

P(a)=A 3 [6 a , 1 +5 a , 3 ]+A 4 [8 <Ti2 + S aA } (a = 1,2,3,4) , (10) 

% this also means that our formalism in its present form only applies to the study of K-edgcs. 
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(fc|V|n) = <n|V|fc) = E (3(a) [ 7 +(a)5 fc , n+ff + 7 ~ (a)S k , n - a ] , (11) 

(7=1 

where we have explicitly used the hermiticity of V, and defined 



>.! 



+ t \ (n + o-)\ , _ I ./ 7 ^rra (n,f) n > a , . 

otherwise. 



The coefficient a ± (n, cr) are reported in table 1. We note that the coefficients 7n 
and (n, a) satisfy to the following relations 

=ln a + (n-cr,a) =ar{n,(j) ^ 

ln+M) =7rt a-{n + a,a) =a+{n,a) . 

Using eq. (11), (12) and (13), we can work out the second-order corrections to the 
energy levels AE n and the normalised perturbed wavefunctions |n) + \n') in the 
usual way (see e.g. ref. [16]). After a somewhat lengthy calculation, we get 

A£;„/^ = 3A 4 (2n 2 + 2n+l) + ^^[( 7 -( ( T)) 2 -( 7 +(a)) 2 ] . (14) 

(7=1 

\n') = E ^ bn (*)l« - *) ~ lM\n + *)] + 

(7=1 

E SrSy W (^£->')i n + ^ + + 7„-(^)7„--<7K)in - (<r + a))] - 

<7,(T' = 1 

E ^fX] bMi^O ^)> + ln(.°)lt-M)\n + 0-))] 



Tiuj a 2 

(7=1 



^|n)E^ [(7„ + W) 2 + (7„"W) 2 ] , (15) 

(7=1 

where AEn" 1 — 3Hlu h±(2n 2 + 2n + 1) arc the first-order corrections to the energy 
levels. 

We are now able to write down explicitly the expression for the QPDF from 
eq. (3). We have 

EKM + <(r)] 2 e-[^ 0) + ^ n \/k B T 

9(r,T) = ^ , (16) 

j- e -[E^ +AE n ]/k B T 

n=0 

where u n (r) = (r\n) are the unperturbed cigenfunctions, i.e. the eigenvectors of the 
one-dimensional harmonic oscillator, and u' n (r) = (r\n') are the corrections (15). 
We have explicitly indicated the truncation of the summations as the integer n M . In 
the computations one has to fix n M by requiring that the corresponding normalised 
Boltzmann factor z nM /Z (Z being the partition function) is negligible up to some 
specified tolerance tol = 10~ M , i.e. 

ZnM =e-^°l+ AE ™]/k B T< tol 
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This condition fixes the number of levels which are included in the perturbative 
series (16). Of course, one also has to check that the energy of the highest level 
included is small compared to some estimate of the the potential well depth V Q . 
If we express energies in eV and temperature in Kelvin, from eq. (17) we get the 
condition 

[M log 10] T x 10~ 4 < V . (18) 

It follows that, for potential well depths of the order of 1 eV, the condition (18) is 
fulfilled for M = 4 — 5 for temperatures up to « 100 K. 



3. Validity and sensitivity of the QPDF 

The above described procedure to build the QPDF relies on two basic hypotheses: 
(i) the classical approximation of lattice vibrations must be inadequate so that 
the quantum treatment of the two-body problem holds, and (ii) the deviations 
of the absorber-neighbour potential from the harmonic approximations must be 
satisfactorily described by a perturbative treatment. As to the validity of condition 
(i), the ratio 1Z Q — fiujjk B T (oj = y/k 2 /fj,) provides a good qualitative indicator: if 
TZ Q is of order one, the quantum energy scale is comparable with the thermal one, 
and the system requires the full quantum description. Regarding condition (ii) , we 
introduce the following parameter 



(19) 

where (. . .)t is the configurational average performed with the QPDF. The indicator 
1Z E gives a measure of the relative strength of the perturbed and unperturbed 
energies, and can always be computed a posteriori in order to assess the validity of 
the perturbative approximation. On the other hand, we also have to be concerned 
with the sensitivity of the QPDF to the parameters of the model potential, i.e. the 
minimum detectable anharmonicity within the present model at fixed temperature 
and potential stiffness (k 2 ). We shall here introduce a simple procedure for assessing 
the QPDF sensitivity, based on the statistical F-test. 

Let us suppose that we are fitting TV experimental data points (ki,xT v ), 
i = 1, . . . , N, to a model that has p adjustable parameters Xj, j = 1, . . . ,p. The 
model predicts a functional relationship between the measured independent and 
dependent variables 

x(fc) = x thc (fc;Ai,...,A p ) . 

In the spirit of the maximum-likelihood method, we want to minimise a fit index 
(or residual function) of the kind 

N 

T = " x th '(h; Ai, ... , A p )] 2 w t (20) 

i=l 

where the w^s are some weight functions. In general, if the standard deviations 
<7j of the experimental data are known independently, the correct choice would be 
Wi cx 1/of. However, depending on the particular algorithm used to perform the 
fit, some other weighting functions may be preferred. 

As a consequence of introducing the third- and fourth-order nonlinearities in 
the model potential, two more floating parameters are available to fit the EXAFS 
spectrum, namely k% and k^. In general, this will cause per se a reduction of 
the best-fit index minimum with respect to the harmonic model. Such a situation 
typically arises in EXAFS data analysis when it is to be decided whether a spectrum 
needs the introduction of an additional shell to be fitted (some physical information 
still present in the data) or, more generally, whether the improvement achievable 



^ d g 2yi(V((r-ro)2)r) = J_ 
k 2 ((r-r ) 2 ) T k 2 



fW<(r-ro)2> r + iM(r-r ) 2 >T 
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by incrementing the number of free parameters is statistically meaningless [17, 18]. 
In our case, we are interested in assessing the sensitivity of the QPDF in capturing 
real physical information regarding the higher-order terms in the Taylor expansion 
of the potential. 

A general theorem in statistics states that the minimum of the residual function 
F min is distributed as a x 2 distribution with v\ — N—(p+l) degrees of freedom [19] 
§. Let us assume we want to compare an harmonic model potential (fit with p 
parameters) to a potential obtained by adding some anharmonicity (fit with p + p' 
parameters, p' = 1,2). We want to assess whether the latter model significantly 
improves the fit (given the automatic improvement following the introduction of 
any additional free parameter) . Let us consider the ratio of the normalised minima 
of the two corresponding fit indeces 

F ~T m up+ P ')H ' (21) 

where = [N — (p + p' + 1)]. It can be shown that, as a consequence of the above 
theorem, the ratio (21) follows an F V1V2 distribution [19], whose density function is 

V^{F)dF= r[(, 1+ , 2 )/2][(, 1 /, 2) F]^ 

In particular, if both the harmonic and anharmonic models were appropriate to 
explaining all of the signal, one would expect the function 



P' \F m in(P+P') 

to follow an F p >^ 2 distribution. The F-test is then conducted as follows. (1) Based 
on some estimate of the experimental standard deviation, one generates an artificial 
experimental data set from the anharmonic model, e.g. by adding some Gaussian 
noise. (2) The two residuals F min {p) and F min (p + p') are calculated and the value 
of / obtained. (3) One can now fix the preferred confidence level c and compare / 
with F£, ^ , given by the following relation 

c=l- [ Fp '-" 2 V p ,, V2 (F)dF . (24) 
Jo 

Here c is the percentage probability of obtaining a reduction F min (p) — F min {p + p') 
as large as actually observed, when the added anharmonicity is not physically 
meaningful. We shall call such confidence level the rejection probability, which 
expresses the probability that the harmonic and anharmonic model are equivalent - 
that is, if c = 1 the two models are by all means statistically indistinguishable, while 
if c = the probability that the harmonic model could explain the data instead of 
the anharmonic one vanishes identically 



4. The Ag— I potential in Silver Iodide 

Silver iodide (Agl) is known for being a highly-anharmonic material [21, 22, 23]. 
Hence, it appears a good candidate for providing a bench-mark to test the above 
described data analysis framework. As an example, we analyse here an Ag K-edge 
EXAFS spectrum collected at T = 77 K. The details of the experiment and of the 

§ This result strictly holds when (i) the measurement errors are normally distributed, and either 
(ii) the model is linear in its parameters or (Hi) the sample size is large enough that the 
uncertainties in the fitted parameters do not extend outside a region in which the model could 
be replaced by a suitable linearised model. Fits in EXAFS are usually at the limit of validity of 
condition (Hi). However, all the statistical analysis that stems from this basic theorem is routinely 
applied in EXAFS data analysis (see e.g. [20]). 
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extraction of the EXAFS signal X exp (k) from the raw absorption data are reported 
elsewhere [25]. 

The classical expression (2) of the PDF has been recently used in an EXAFS 
study of Agl at the I K-edge to measure the first three coefficients of the Taylor 
expansion of the Ag-I potential at T = 300 K and T = 600 K [24]. Following ref [24], 
we get TZ Q (T = 77) « 0.5. We conclude that at T = 77 K the classical expression of 
the PDF is non longer valid and one has to work in the quantum regime. In order 
to calculate the model EXAFS signal we re-write eq. (1) by using the standard 
formula of single shell single-scattering EXAFS in the following fashion 

r+oo -2[r +u]/\(k) , . 

X (k, r) = S 2 o m J g(u, T) - Im {/i(7r, k) e 2 ^ e 2 ^o+«] j du , (25) 

where g(u,T) is the QPDF as calculated from cq. (16) and we have made the 
substitution u = r — r . Here Nj is the coordination number of the I ions, which we 
fixed at the crystallographic value N\ — 4, <5A g is the central atom phase shift and 
fi(ir,k) = \fi(iT, k)\ exp(i^i) is the complex backscattering amplitude of the I ions. 
The constant S% is the usual reduction factor which accounts for the inelastic losses, 
which we fixed at the value S% — 0.73 [25]. The photoelectron mean free path is 
here assumed to depend on the wavevector k as A = k/t] [2], where rj is a constant. 
The backscattering amplitude fi(ir,k) and total phase shift Aip = 25 \ g + <pi have 
been taken from the tables in ref. [2] (reproduced from calculations based on the 
Herman Skillman wavefunctions [26]). 

Following the arguments developed in section 3, we carry out two separate 
non-linear least-square fittings, by comparing the experimental EXAFS X exp (k) 
with the model signal calculated by eq. (25) from the two separate sets of floating 
parameters {A} = {ro, AEq, rj, k^\ and {A'} = {ro, AEq, r], k2, k%, k^}. The floating 
parameter AEq must be included in the fit as usual in order to compensate for the 
uncertainty associated with the true value of the threshold energy, i.e. the minimum 
energy required to free the photoelectron. 

The best-fit values of the free parameters are reported in table 2 for both 
the harmonic and anharmonic models, alongside with the corresponding fit index 
minimum J- m in- The quality of the fit with the anharmonic model is shown in 
fig. 1 (a), while the corresponding effective Ag-I potential is drawn in fig. 1 (b). In 
table 3 we report for comparison the values of the potential parameters as measured 
in ref. [24]. We see that the overall agreement is good, although the value of fc 4 
reported in ref. [24] does not seem to be reliable. Moreover, if we use eq. (19) to 
estimate the strength of the perturbation energy corresponding to the best-fit values 
of kz and k&, we get 1Z E w 0.12. We conclude that our results of the analysis with the 
QPDF are consistent with the perturbative hypothesis. It is instructive to observe 
that the same analysis performed using the classical expression of the PDF eq. (2) 
with the parameter set {A'} always yields a vanishing ki at the minimum T m in- 
Correspondingly, the fourth-order constant k^ is raised to unphysically high values, 
independently of the initial guess of the set {A'}. This scenario clearly confirms the 
inadequacy of the classical treatment of the pair dynamics in the present case. 

We turn now to assessing the validity of our treatment on statistical grounds. 
The value of / corresponding to the reduction of T following the introduction 
of &3 and &4 can be calculated by eq. (23). Substituting N = 150 and p' = 2, 
we get / w 16.6. By substituting in turn F p c , v ^ = f in eq. (24) we get the 
corresponding rejection probability c«3x 10~ 7 . We arc then allowed to conclude 
that the anharmonic model is here capturing a real physical feature of the Ag-I 
pair dynamics. It is instructive to demonstrate this conclusion in a more pictorial 
fashion. It is a simple corollary to the first theorem mentioned in section 3 that, if 
only v parameters are varied while keeping the other p — v fixed at their best-fit 
value, the function AT = T{\\, A2, . . . , X v ) — T m i n follows a \ 2 distribution with 
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v degrees of freedom. When v — 2 this result provides the errors on selected 
parameter pairs (Ai,A2) in the form of confidence ellipses through the simple 
condition AjF(A 1; A 2 ) = [x 2 ]£=2- Here [x 2 ]£=2 is the value of the x 2 variable 
corresponding to the required confidence level (i.e. rejection probability) c for v = 2. 
In fig. 2 we show the contour levels of the function AT = T(ki, k^)—T m in computed 
for the minimum obtained within the harmonic model. The presence of a significant 
region (corresponding to the 99 % confidence level [x 2 ]°=2 = 9.21) of negative values 
away from the computed minimum explains the dramatic improvement of the fit 
upon introducing the nonlinearities in the potential. The same confidence analysis 
for the anharmonic model performed in four different parameter subspaces {Ai, A 2 } 
is reported in fig. 3, showing the quality of the best-fit minimum. 

We end this section by showing how one can conduct the F test described 
in section 3 to examine the QPDF sensitivity in the present case. Let us fix the 
temperature and the harmonic constant fc 2 at its best-fit value. We can then use 
eqs. (25), (20) and (23) with p' = 1 to calculate the value of / for any choice of 
k 3 and fct, where in place of the experimental spectrum we use an artificial data 
set constructed as described in section 3. In particular, we just add to the model 
signal a Gaussian noise with standard deviation a{k,T) = cr [ko/k] 3 ' 2 \JT /To, with 
al = 0.0016, fc = 12.7 A" 1 and T = 77 K (see fig. 1 (a)). Finally, we calculate 
the corresponding rejection probability by means of eq. (24) with F£, ^ = /. For 
example, this procedure can be used to construct the functions 0(^3,^4 = 0,T) 
and c(fc 3 = 0, &4, T) (one-parameter sensitivity curves). Alternatively, the same 
procedure with p' = 2 can be used to look at the contour sections of the function 
c(ks, &4, T) (two-parameter sensitivity curves). In fig. 4 we show an example of one- 
parameter sensitivity curves calculated at three different temperatures for both the 
&3 and &4 parameters. We clearly see that the lowest detectable anharmonicity is 
well below the measured one in both cases. In particular, the rejection probability 
decays exponentially with increasing magnitude of the anharmonic constant. 

5. Conclusions 

In this paper we introduced a pair distribution function valid in the quantum 
regime based on a Taylor expansion of the absorber-neighbour potential, suitable 
for the analysis of single-scattering EXAFS. In particular, we used ordinary time- 
independent non-degenerate quantum perturbation theory to cast the QPDF in a 
simple analytical form, which can be easily calculated in a subroutine incorporated 
in the fitting program. Moreover, we showed how the limits of validity (sensibility 
of the quantum treatment and perturbative hypothesis) can be monitored and we 
described how to implement a simple statistical test to estimate the sensitivity of the 
QPDF to the third- and fourth-order terms in the potential. The latter procedure 
can be used a posteriori to check whether the best-fit values of the anharmonic 
parameters are above the sensitivity level (minimum detectable anharmonicity). 
Alternatively, the same procedure may be applied a priori in order to assess whether 
the QPDF is suitable for the analysis of the problem at hand. We applied our 
formalism to the case of Silver Iodide, showing how the potential anharmonicity 
can be measured in a temperature range where the classical expression of the PDF 
could not be used. In particular, we demonstrated how a simple harmonic model 
is not adequate to describing the dynamics of the Ag-I pair, in agreement with 
previous studies performed in the classical regime. 

As a final remark, it should be noted that the example we chose to test the 
QPDF concerns a very symmetric lattice structure. It is well known that an 
additional broadening of the distributions of absorber-neighbour distances may 
be produced by static disorder. The latter may be associated for example with 
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the presence of a coordination shell made of N identical atoms at slightly different 
distances from the photoabsorber, that can not be resolved in different subshclls 
(e.g. Ni at distance r\ and N 2 at distance r 2 , with Ni + N 2 = N). For small static 
disorder (|n — r 2 \ <n,r 2 ), one can prove that such shell is equivalent to a shell with 
coordination N and mean distance ro = (N\ri + N 2 r 2 )/N, provided one introduces 
in the Debye- Waller factor the additional variance <7g tat = N\N 2 \r\ — r 2 \ 2 /N 2 [2]. In 
the framework of our model, this is expected to correspondingly rescale the harmonic 
constant k 2 . However, it should not alter in general the information carried by the 
EXAFS signal regarding the anharmonic terms of the absorber-neighbour effective 
potential. Hence, we expect that our model of lattice anharmonicities may be 
sensibly used in its present form also in the presence of small static disorder. 

Concluding, we developed an easy technique to study interatomic potentials 
from single-scattering EXAFS in the quantum regime of atomic vibrations. A full- 
featured FORTRAN package containing the relevant programs is made available by 
the author to the interested scientists upon request. 
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Table 2. Best-fit values of of the fit index (20) and parameters describing 
the Ag-I potential. The weight functions used are here u>i = 1/^/50 Vi. 



model 


T 

^ mm 


fc 2 (eV A- 2 ) 


k 3 (cV A- a ) 


fc 4 (cV A- 4 ) 


n (A- 2 ) 


(A) 


AE (cV) 


{A} 


41.0 
31.5 


2.11(3) 
1.86(3) 


-7.6(4) 


49(6) 


0.84(5) 
0.82(5) 


2.86(1) 
2.87(1) 


-48(1) 
-43(1) 



Table 3. Ag-I potential parameters as measured in ref. [24]. 



T(K) 


fc 2 (cV A- 2 ) 


fc 3 (eV A- 3 ) 


fc 4 (eV A" 4 ) 


300 


2.4(1) 


-5.0(1) 


0.2(0.7) 


600 


2.7(1) 


-5.3(1) 


0.08(0.4) 




(a) (b) 

Figure 1. (a) Upper frame: Ag K-edge EXAFS signal at T = 77 K (symbols) 
and fit obtained with eq. (25) and parameter set {A'}. Lower frame: squared 
residuals (symbols) and fit with the law u 2 (k) = <Tg [ko/k] 3 (solid line), with 
a\ = 0.0016 and k = 12.7 A -1 , (b). Effective Ag-I potential. Also shown 
are the harmonic levels used to calculate the QPDF, and the corresponding 
normalised Boltzmann factors z n /Z. 
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Figure 2. Contour levels of the function AjF = T(k2,ki) — T m i n for the 
minimum corresponding to the harmonic model. The contours correspond, 
from left to right, to A^=-[ X 2 ]™ 2 , -[x 2 ]°^ 2 ,-[x 2 ]°^ 2 , [xX=2>ixX=2> 
[X 2 ]°i 2 and [ X X °l 
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Figure 3. Contour levels of the function AJF = JF(Ai,A 2 ) — T m i n for 
the minimum corresponding to the anharmonc model. Four combinations 
of floating parameters pairs {Ai,A 2 } are shown. Contours correspond to 
confidence levels 50%, 70%, 90%, (95%) and 99 %. 
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Figure 4. Rejection probability as a function of the potential anharmonicity 
calculated by eq. (24) at different temperatures for fixed harmonic constant. 
Parameters are N = 150, \x = 9.74 X 10 -26 Kg and k 2 = 1.86 eV/A 2 . (a) cubic 
nonlinearity (b) quartic nonlinearity 



